Close contact fluctuations: the seeding of signalling domains in the immunological 

synapse 
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We analyse the size and density of thermally induced regions of close contact in celbcell contact 
interfaces within a harmonic potential approximation, estimating these regions to be below 1/10 th 
of a micron across. Our calculations indicate that as the distance between the close contact 
threshold depth and the mean membrane-membrane separation increases, the density of close 
contact patches decreases exponentially while there is only a minimal variation in their mean size. 
The technique developed can be used to calculate the probability of first crossing in reflection 
symmetry violating systems. 
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Surface contact between cells is a key mechanism for in- 
formation transfer in many biological systems. These 
can be both long term or permanent connections as in 
the neurological synapse, or as discovered more recently, 
transient and highly dynamic as in the immunological 
synapse [l| . T-cells (a class of lymphocytes) make tran- 
sient contact with 'target' cells whilst scanning for the 
presence of their specific antigen, antigen recognition re- 
sulting in the stabilisation of the contact and genera- 
tion of a macroscopic receptor patternation in the con- 
tact interface, or a so called immunological synapse 
A fundamental observation is that the contact interface 
is heterogeneous, both in the physical separation of the 
two cell surfaces Q and in the local signalling proper- 
ties P, [3, 0] . Differences in the extracellular lengths of 
key molecules/bonds is believed to underpin both these 
processes with a predominant division between short and 
long bond length molecular species. Of note is that es- 
sential antigen signalling receptors, such as the T-cell re- 
ceptor (TCR) , are small molecules with a ligand-receptor 
bond length of 14nm (membrane to membrane span) [l[ , 
while an essential phosphatase (CD45), a major compo- 
nent of the glycocalyx, has a length of 25-40nm and is 
not believed to have a natural ligand. T cell signalling, 
or antigen detection, thus requires tight cell:cell contact 
to allow TCR binding, whilst such regions necessarily 
require the spatial exclusion of the large molecules com- 
prising the glycocalyx. Spatial heterogeneity in the mem- 
brane profile within the contact interface is therefore es- 
sential for the functioning of the cell contact. Early pat- 
terns (50sec) in cell interfaces show random small clusters 
of TCRs [3, H|, regions where signalling intermediaries 
appear to congregate. These regions of close contact are 
presumably formed from fluctuations in the initial con- 
tact surfaces. At later times signalling appears to be 
focused in distinct stable microclusters [5|- This depen- 
dence of signalling on spatial heterogeneity introduces a 
key 'exposure' problem; ligand detection requires that re- 
gions of close contact comprise a significant area within 
the interface while they must be sufficiently large that 
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FIG. 1: A ID profile of a 2D membrane fluctuating around 
the line < Z >= 0, from simulated data. X+ & X- are the 
respective sizes of the patches above and below the threshold 
depth Z = A = lOmn. 



they can be stabilised when segregation is energetically 
favourable Q . We examine the spatial statistics of these 
regions of close contact using a linear stochastic model 
for thermal fluctuations of the membrane separation. 

In this letter, our interest is in the size and den- 
sity of regions of close contact (eg membrane-membrane 
separation <20 nm) where effective TCR ligand bind- 
ing can occur. We utilise a linearised version of the 
synapse reaction-diffusion equations to model pre 

pattern dynamics, reducing to a single equation for the 
membrane-membrane separation Z around a mean sepa- 
ration (25-50nm) determined by the glycocalyx potential 
and receptor-ligand bond equilibrium. In this regime, the 
fluctuation Z{x, t) has dynamics 
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where B is the membrane rigidity, r the surface tension, 
M the membrane damping constant and A parametrises 
the rate of relaxation of the membranes close to equi- 
librium, i.e. the strength of the harmonic approx- 
imation to the potential well. The thermal noise 
rj(x,t) is defined using a fluctuation-dissipation relation 
< r)(x,t)ri(x',t') >= 2k B TM5 2 {x-x')5{t-t'), x,x' be- 
ing points in the contact interface. The solution Z(x, i) is 
a Gaussian variate. We wish to calculate the probability 
that the displacement Z(x, t) lies below a 'close contact' 
threshold —A where A ~ 5 — 30 nm is the membrane- 
membrane displacement from the mean required for ef- 
ficient TCR binding. We identify the region Z < — A 
as a region of close contact and determine the average 
size of these close contact patches. A point to note is 
the symmetry violation of the system around Z = —A; 
specifically the average size of a patch above this line 
(designated by +) is not the same as one below this line 
(designated by — ). The statistics for Z < — A are identi- 
cal to those for Z > A; thus for presentation we will use 
Z > A as the threshold. 

We start by defining the sign (conditional) correla- 
tor for an arbitrary displacement x in the contact in- 
terface (relative to the origin) [j| Q A + =< sgn[Z(x) — 
A ] >z(6)>A and,4_ =< sga[Z(x)-A] > z( g )<A , < .. > F 
denoting the average over states where condition F holds. 
For simplicity we assume Z is in stationary equilibrium 
and thus initial conditions can be ignored. Z(0), Z(x) 
define a two variable joint Gaussian probability distribu- 
tion with zero means, variance c\\ =< Z 2 (0) > and co- 
variance 012(2?) =< Z(0)Z(x) >. By translational sym- 
metry the covariance matrix and A± only depend on the 
spatial displacement x — \x\ between the membranes. 
Thus we drop explicit reference to x for simplicity. The 
symmetry relation A + (x, A) = — A—(x, — A) means that 
only A + needs to be evaluated. 

An ensemble averaging over the two- variable Gaussian 
distribution gives 
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where the lower limit follows from the condition Z(0) > 



A. Here detc = 



c?o and the normalisation con- 



stant N + is defined by the error function 
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which is in fact the probability of ob- 



serving a separation Z(x) > A at an arbitrary point 
x. We define the patch sizes X± for regions where 
Z > A, Z < A respectively (in 2D along an arbitrary 
vector), Fig. [T] To evaluate < X± >, we need to 
evaluate A' ± (0, A) where the prime refers to a deriva- 
tive with respect to 1 = |a?|. This follows from the 
relation A±(x, A) = 1 — 2x/ < X± > as separation 



x — > 0, a consequence of the fact that the probability 
of finding a crossing (i.e. Z — A) in a small interval 
of length x is <j ^ E ±> ■ This gives us the exact relation 

< x ± >= ~3£{oy 
We proceed to compute the derivative as A'±(x) = 

Itf-^ Eq. © giving 
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The relevant correlator in 2D is given by 
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x = ( 4 x and Kq is a modified Bessel function of de- 
gree 0. The final integral uses a Bessel function identity 
[lCl ]. Therefore cn 



and for small x wc find 
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A + (x, A) ~ 1 — i(log e x)?C, C a constant. Thus cross- 
ings fail to conform to the assumptions above, specifi- 
cally C12 is not twice differentiable at x = This 
is a familiar consequence of Brownian motion crossing 
behaviour and stems from the high frequency noise com- 
ponent of 77 that causes repeated crossing of the threshold 
in between large excursions away from the threshold. We 
regularise the divergence by introducing an infra-red cut- 
off in the noise, thus correlator ((4]) becomes 
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where Jo is a Bessel function of degree and cut-off k m — 
2tt / e is given by the smallest length scale e in the system. 
This length scale is on a sub nanometer scale, e.g. the 
width of lipid molecule head in the membrane. A regular 
expansion for cyi at small x now follows 



Provided e is sufficiently small we have a consistent 
regularisation with c" 2 < at x — 0. We thus obtain the 
following expression 

M A 1/4 

A' + (0,A) = -2-± (-) exp(-A 2 /2 Cll ) (7) 



where / = 47ry cf>/ [log sinh0^ depends only on 

system parameters. We have retained only the leading 
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FIG. 2: Variation of < X+ > against A for Z > A: theoret- 
ical estimate from Eq. ((8} vs numerical simulation taken on 
a lattice, size 1000, spacing lnm, over 1024 runs. A coloured 
noise spectrum is used, derived by projection from 2D, giving 

< r)(x,t)rj(x',t') >= 2k B TMs{x - x')6 2 (x - x )S(t - t') with 

S(k) = (M^hlf 4 JL where £ (fc) = l og [ r+2B^ _ 

W 4coshi0(fc)' ^ v ' aL ^/4BMa{k) 

l ( T +2Bk^) 2 7 
V 4,BMa(k) J' 

order in the cut-off for simplicity. The mean sizes of the 
patches above and below the line Z = A now follow, 

< X± >= ^(|) 1/4 exp(A 2 /2 Cll ) (8) 

where the normalisation constant iV_ is defined as 
N-^ 1 = 1— JV+ . The dependence on the cut-off is weak 

while the length scale is determined by ^JcuB/ksT. 
For a symmetry preserving system with A = we have 

< X± > = |-(-r) 1//4 - Suitable values for the system pa- 
rameters are p : B = 11.8 ksT, r = 5650fcsT ^m~ 2 , 
M = 4.7 x 10 6 k B T s ^m~ 4 and e = lnm, while 
A = 6.0 x 10 ksT [ivnT^ is approximated from the lin- 
earised reaction-diffusion equation as rx CD45 density, 
the latter being approximately 100 molecules /im~ 2 . This 
follows from the force expression in synapse reaction dif- 
fusion equations, J~J. k(z — li)Ci, a sum over all molecules 
Ci that impose a force on the membrane (bond length lj) 
with a spring constant k ~ r, [l,[l3]- In early signalling, 
CD45 will be the dominant component. These values 
give / = 2.5, < Xfc(A) > |a=o ~ 84 nm. The variation 
of < X± (A) > with A is illustrated in Fig. \2b Density 
fluctuations in the Ci will causes fluctuations in A which 
can be included as a "non-equilibrium temperature" in 
Eqn. (1) (fluctuation-dissipation relation); however this 
is beyond the current minimalist model. 

As the threshold A increases above zero the regions 



FIG. 3: Probability density of humps around the line Z — A 
as a function of the average threshold A (nm) as predicted 
from Eq. © 



Z > A develop into isolated patches in 2D. We can use 
the mean size to estimate the patch density phumps by 
a mean field approximation phum P s< X + > 2 — N^ 1 to 
obtain 

, N+ . A N 1/2 . A 2 

Phumps = -JJ- (-g) exp(- — ). (9) 

The expected decline in the density of patches as A in- 
creases is shown in Fig. [3l For large A 3> en the 

i 
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leading behaviour is < X + (A) >^ / , ^ A (^) and 

a / \ A 2 
Phumps ~ f J exp ( - £^ ) . These asymptotic ap- 
proximations capture the contrasting weak decline of the 
width < X + >, and strong decay of the hump density 
Phumps with A in Figs. In ID, regions with Z > A 

are always disconnected so the patch density Phumps can 
be defined for all values of A. Further, in ID, there are no 
divergences, whilst in higher dimension the divergences 
are more severe. These properties result from the inter- 
play between the 4 th order PDE Eq. Q] and the volume 
of phase space. 

The probability density function of the distance be- 
tween crossings can also be approximated. This utilises 
the probability distribution for the number of crossings 
of the line Z = A, which is computed by generalising the 
traditional 'persistence' analysis [1, We need to dis- 
criminate between the two types of crossings, a crossing 
from Z>AtoZ<A&sx increases, and the con- 
verse. Let Pn(x) denote the probability that an interval 
of length x contains n crossings of Z(x, t) across the refer- 
ence level A when Z > A at the extreme left, and p~ (x) 
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is the corresponding probability with Z < A at the ex- 
treme left; in 2D we consider moving a distance x along 
a specified vector. Then under an independent interval 
approximation, for n > 1, their Laplace transforms have 
the forms 
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for n odd and even respectively. Here P± {x) is the prob- 
ability density for the distance x between crossings (+ : 
Z > A, - : Z < A), and < X > = (< X+ > + < X^ > 
)/2 is the average distance between consecutive crossings 
(any type). Using the identities Y^=v Pn(X) = 1, we 
can now show that On = s^ 1 -N±(l~P±)/(2s 2 < X >), 
which agrees with [9( when A = 0. Employing the iden- 
tity A±{x) — X^^Lo ( — 1)™ P±( x )i we arrive at two cou- 
pled equations relating P± and A±, 
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Solving these equations then gives the desired pdfs 
P ± (x). 

To summmarize, for the harmonic potential membrane 
model we have an exact analytic calculation for the mean 
size of close contact patches, < X + (A) >, our calcula- 
tions suggesting that these are on the scale of tens of 
ran. The scale is primarily determined by the combi- 
nation 



I cuB / ksT and has a leading order behaviour 

going as l/A for large A. This small patch size im- 
plies that multiple receptor bindings are unlikely within 
a patch and close contact patches are unobservable by 
traditional light microscopy. The small size also implies 
that phosphatase exclusion (CD45) probably results from 
density fluctuations, ie a specific exclusion mechanism is 
not required in contrast to that needed at larger sizes 
0, 0|- The density of patches decays rapidly with the 
threshold A on a length scale of ^Jc\\ ~ 5.4nm, Eq. (8), 
and indicates that cell membranes must be highly flexible 
otherwise the glycocalyx would impose too large a barrier 
to allow formation of close contact regions (A increasing 
with membrane elasticity). In particular, the glycoca- 
lyx cannot be too deep relative to the size of the TCR 
ligand-receptor bond length (14nm) since otherwise the 



density of patches becomes too small for antigen detec- 
tion. The probability of T cell signalling depends on the 
ability of the TCR to bind it's ligand and is thus crucially 
dependent on the area of close contact regions within 
the cell: cell interface which varies as NT , and on the 
size of those close contact patches. There is an enhance- 
ment in triggering as patch sizes increase above 150nm 
fl3j |: thus our estimates suggest that early signalling re- 
lies on patches below this size and enhancement effects 
only occur upon aggregation and stabilisation of clusters 
as the immunological synapse forms. Such conclusions 
are somewhat reminiscent of [l2[ where formation of a 
synapse was related to a critical value of the system pa- 
rameters (albeit without evaluating the patch size). We, 
however, go beyond such qualitative predictions. Our cal- 
culations clearly suggest that the membrane correlation 
length is a determining factor in the area of close contact 
regions in the interface, which with our parameters limits 
the threshold to A < 16 nm. 
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